Far Infrared Properties of M Dwarfs 



Thomas. N. Gautier III ^, G. H. Rieke ^, John Stansberry ^, Geoffrey C. Bryden ^, Karl R. 
Stapelfeldt ^, Michael W. Werner ^, Charles A. Beichman ^, Christine Chen ^, Kate Su ^, 
David Trilling ^, Brian M. Patten ^ and Thomas L. Roellig ^ 

ABSTRACT 

We report the mid- and far-infrared properties of nearby M dwarfs. 
Spitzer/MIPS measurements were obtained for a sample of 62 stars at 24 fim, 
with subsamples of 41 and 20 stars observed at 70 /im and 160 /^m respectively. 
We compare the results with current models of M star photospheres and look for 
indications of circumstellar dust in the form of significant deviations of K-[24 fim] 
colors and 70 /im / 24 yum flux ratios from the average M star values. At 24 fim, 
all 62 of the targets were detected; 70 fim detections were achieved for 20 targets 
in the subsample observed; and no detections were seen in the 160 fim subsam- 
ple. No clear far-infrared excesses were detected in our sample. The average far 
infrared excess relative to the photospheric emission of the M stars is at least four 
times smaller than the similar average for a sample of solar-type stars. However, 
this limit allows the average fractional infrared luminosity in the M-star sample 
to be similar to that for more massive stars. We have also set low limits (10~^ 
to 10~^ MEarth depending on location) for the maximum mass of dust possible 
around our stars. 

Subject headings: Infrared: stars - stars: late-type - stars: fundamental parame- 
ters (colors) - planetary systems 
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Introduction 



The study of dust around main sequence stars addresses intriguing questions in basic 
astrophysics while providing insights into the possible existence of planetary systems around 
these stars. The pion eering results from IRAS on infrared excess around stars like Vega 
( lAumann. et al.lll984j ) rapidly led to the realization that the circumstellar dust responsible 
for this emission could not survive for the lifetime of the star and that it must be replenished 
- from cometary or asteroid-type objects - just as the zodiacal cloud in our Solar System is 
replenished. The maintenance of this dust ultimately depends on collisions between asteroid- 
size objects, yielding fragments that collide with each other and eventually grind down to fine 
dust that is heated by the star and appears as an infrared excess. Thus the characteristics 
of circumstellar dust around a main sequence star (referred to hereafter as a "debris disk") 
reflect the planet-building and planet-destroying processes occurring around that star. 



The Spitzer Space Telescope ( IWerner. et al.ll2004l ). building on the foundation estab- 
lished by IRAS and ISO, has already produced major advances in our understanding of 
the debris disk phe nomenon. These include demonstration of the stochastic nature of de- 
bris disk evolution (IRieke. et al.ll2005l : ISu. et al.ll2006l : iBeichman. et al.ll2005al ). the first de- 
tection of debris disks around F, G, and K stars with known radial velocity companions 



( iBeichman. et al.l l2005b[) and investigations of the frequency of occurence of debris disks 
around solar type stars (IKim. et al.ll2005l : iBryden. et al.ll2006l ). Here we present the results 
from a Spitzer search for dust associated with M dwarfs. 

M dwarfs, with masses between 0.08 and 0.5 Mq, are the most common stars in the 
galactic plane. Thus the study of debris disks around them is needed to complete our 
understanding of the nature and incidence of debris disks in the solar neighborhood and 
the dependence of planet-forming processes on the mass of the central star. Since M stars 
are the most numerous stellar type, measurement of their debris disk frequency has large 
implications for the overall population of planetary systems in the galaxy. It may be that 
the nearest extrasolar planets to the Sun are orbiting M stars, a possibility that would be 
supported if M dwarfs frequently host debris systems. 

However, the timescales and processes for du st expulsion may be co nsiderably different 
for M stars than for their more luminous cousins (jPlavchan. et al.ll2005l ). Also, early results 
from precision radial velocity surv eys have suggested that M stars have a lower frequency of 
giant planets than F and G stars (iButler. et al.l 120041 ). It is important to determine if such 
a marked difference is also present in the relative frequency of debris disks between the M 
dwarfs and earlier types. 



Finally, besides extending our understanding of the evolution of planetary systems by 
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investigating the frequency of debris disks around M dwarfs, these far-infrared measurements 
can determine basic parameters of the M stars themselves. We compare our data with current 
models of M dwarf photospheric emission and demonstrate how to refine determinations of 
M dwarf effective temperatures with the accurate 24 /im data available with Spitzer. 



2. New Observations 

We observed 62, 41, and 20 M dwarfs, respec tively, at 24 , 70, and 160/im, using the 



Multiband Imaging Photometer for Spitzer (MIPS; iRieke. et al.ll2004j ). Our target stars are 
listed in Table [H Stars of spectral type M7 and earlier (a total of 40), plus one M9 star, are 
part of a MIPS program to obtain photometry of all stars lying within 5 pc of the sun {Spitzer 
program ID 52). Twenty one later- type stars (M6 - M9) lying at distances of 6 to 18 pc were 
generally selected to be the brightest members of their spectral classes {Spitzer program ID 
56). These later- type stars were observed only at 24 fim. While age determination of field 
M stars is notoriously difficult it is likely that the stars in our sample are nearly all older 
than 1 Gy. 

The 24 fim observations were made using the MIPS photometry observation template 
with the small dither size. For the nearby stars we generally used 1 cycle of the dither 
pattern (48s integration time). IL Aqr, LHS 3685 and LP 944-20 were exceptions receiving 
180s, 136s and 1500s integration time respectively. This strategy typically resulted in very 
high signal-to-noise ratios (SNRs) at 24 /im. The more distant, later type, stars received 
from 1 to 9 dither cycles (48s to 1380s) at 24 yum; many of them were measured only to SNR 
~10. 

The 70 /xm observations were also made using the MI PS photometry ob servation tem- 
plate with the small dither size and the default image scale (IRieke. et al.ll2004j ). In this band, 
the required integration times were much longer than at 24 fim and therefore automatically 
resulted in several cycles of the dither pattern. We limited the maximum integration time 
at 70 /xm to ~ 1000 seconds, both to fit the program within its time allocation and be- 
cause the noise behavior of MIPS at longer integration times was unknown at the time these 
observations were scheduled (before the launch of Spitzer) . Our goal of SNR ~ 6 on the pho- 
tosphere was reached for many of the 70 fim observations but the limitation on integration 
time resulted in low SNR detections or upper limits at 70 fim for the fainter members of the 
sample. 



To explore the possibility that these low-luminosity stars might possess very cold debris 
disks, 160 fim photometry was carried out for 20 stars in the sample. Detection of the stellar 
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photospheres, whose median flux densities are expected to be a few mJy, was not acfii evable 
given Spitzer's measured extragalactic confusion limit of 8 mJy, la (IDole. et al.ll2004l ). The 
observations were therefore designed only to detect a source 3-5 times the confusion noise 
level. Four cycles of 160 fim photometry, totaling 84 sec of exposure time for each source, 
were used. 

The data were process ed using the M I PS in strument team Data Analysis Tool (DAT, 
version 2.73) described by iGordon. et al.l (120051 ). For the 24 /xm data, basic processing 
included slope fitting, fiat-fielding, and corrections for droop and "jailbar" artifacts (DC 
offsets between adjacent de tector columns that appe ar in measurements of bright sources). 
Further details are given by lEngelbracht. et al.l (120071 ). The 70 /xm data processing also used 
the DAT, but was si milar to that of the Spitzer pipeline version S12. Details are given in 
Gordon, et al.l (120041 ). Mosaics were constructed using pixels interpolated to about 1/2 the 
native pixel scale of the arrays. 

Fluxes were measured using fairly small apertures to improve the signal-to-noise ratios 
of the measurements, particularly given the relatively high spatial variability of the sky 
background at 70 /im. Most targets were measured using apertures 9.96" and 39.4" in 
diameter at 24 and 70 fim. The observed PSF full width at half maximum is 6.4" and 19.3" 
, at 24 and 70 fim respectively. In a few cases nearby sources contributed significant flux at 
the target location, so we used apertures 25% smaller than those just described to reduce 
contamination. In one case (LHS 58 & 59, identified in the notes to Tabled]), we deblended 
two sources at 70 /im using the model PSFs. The ~ 10" aperture at 24 fim corresponds 
to disk radii of 7-90 AU at the various distances of our stars. At 70 /im, where only stars 
within 5 pc were observed, the ~ 40" aperture corresponds to disk radii of 36-200 AU. 

After extraction, the measured fluxes were correc ted to infin ite aperture. Aperture 
corrections were computed using STinyTim model PSFs (jKristll2005l ) for a 3000 K blackbody 
source. The model PSFs were smoothed until they provided good agreement with observed 
stellar PSFs out through the first Airy maximum. 

Targets were usually identified in the 24 /im images using their J2000 coordinates and 
proper motions from SIMBAD. In most instances the agreement between the predicted and 
observed position was excellent (see Notes to Table [1] for cases where there was a significant 
disagreement). However, many of our later M type stars are rather dim, even at 24 /xm, 
and special care was needed to ensure that the correct object was measured among many 
others in these fields. We confirmed a number of identifications with IRAC 3.6 nm images 



from p rogram 35 (Multiplicity and Infrared Colors of Nearby MLT Dwarfs, iPatten. et al. 



(120061 )). The 24 /im photometric aperture was centered on the photocenter of the target. 
We used the same procedure at 70 /im when the source was clearly detected and there were 
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no confusing sources nearby. In cases where the 70 /im centroid position disagreed by more 
than 5" with the 24 fim position, we forced the 70 fim aperture to be centered at the 24 /im 
position. 

The fluxes we report are based on conversion factors of 1.048 yuJy /arcsec^/(DN/s) and 
16.5 mJy/arcsec^/U70 at 24 and 70 /im, respectively (U70 = MIPS 70 /im unit). Table [1] 
lists the observed fluxes, their formal la uncertainties (not including systematic calibration 
uncertainties) and 3a upper limits where no flux could be determined. 

For most of the observed stars, the error in the 24 fim photometry is dominated by cali- 
bration uncertainty, wh ich contributes systematic 24 fim errors at the level of ~5% (see e.g. 
Beichman. et al.ll2005al ). However, we do not depend on the absolute calibration of the mea- 
surements to identify debris disks. We are instead hmite d by the scatter in the predictions 
of photospheric flux density from shorter wavelengths (see lSu. et al.ll2006l ). Background con- 
tamination by galactic cirrus or by extragalactic confusion is of relatively little importance. 
Dim sources, however, are limited by photon noise and by flat fielding uncertainties created 
by latent dark images. The effective noise level for the dimmest sources is ~0.05 mJy. Except 
for the faintest sources, we can detect excesses down to ~15% above the photosphere in this 
band. 

At 70 fim the flux error is a combination of calibration uncertainty and background 
noise. The variable sky background, a combination of galactic cirrus and extragalactic con- 
fusion, presents a noise floor that cannot be improved with increases in integration time. 
This confusion noise is reached with integra tions of lOO's to lOO O's of seconds depending on 
the local level of galactic cirrus brightness (IBryden. et al.ll2006l ). This variation of the sky 
background within each frame, when convolved with our chosen beam size, is listed as the 
70/im flux error in Table 1. On top of this background noise, we have included a 10% system- 
atic error at 70 pm, consistent with the uncertainty found wi thin photospheric measurements 
of FGK stars (IBeichman. et al.ll2005al : iBryden. et al.ll2006l ). A thor ough discussion of th e 
uncertainties in the 70 fim data can be found in these papers and in lGordon. et al.l (120041 ). 
Excesses of 30% above the predicted photosphere for bright stars at 70 /im should be readily 
detectable. 

At 160 /xm, the DAT was used to calibrate the data, and mosaic the images at a pixel 
scale of 8"/pixel (half the native scale). The data were flux-calibrated using the standard 
conversion factor of 42.6 MJy/sr/MIPS160, where MIPS160 are the instrumental units that 
result after the dat a are normalized by th e stim flashes. The uncertainty on the absolute 
calibration is 12% (IStansberry. et al.l 120071 ). The MIPS 160 fim array suffers from a near- 
IR spectral leak that produces a strong ghost image overlapping the true 160 /xm image for 
sources with stellar temperatures and J magnitudes brighter than 5.5 (MIPS data handbook. 
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http://ssc.spitzer.caltech.edu/mips/dh/ ). Only four of our brightest stars have a possibility 
of being effected by this leak. One of these stars (LHS 70) is the single case where one of 
our targets was detected at 160 fim, the spectral leak was detected and not the true 160 /xm 
emission. This detection is ignored and we give 3a upper limits on the 160 fim fluxes for all 
observed targets. 

In most cases these limits are estimated from the scatter of pixels falling within a sky 
annulus with radii of 64" - 128", and normalized to an aperture radius of 24". In some cases 
we also estimated the sky noise by dropping multiple 24" apertures on the images (away 
from the source position), and using the scatter of the fluxes in those apertures instead. The 
limits are corrected for t he finite ape rture size using a correcti on factor of 2.37, derived from 
STiny Tim model PSFs JKristlboOsh fitted to observed PSFs Jstansberrv. et al.ll2007h . The 
resulting 3a upper limits are presented in Table [H There is considerable variation in the 
sensitivity achieved, due to the effects of galactic cirrus emission. The median upper limit 
is 74 mJy, but the limits for individual stars range from 53 to 963 mJy. 



3. Interpretation 

Because the intrinsic mid- and far-infrared properties of M stars are not well known, we 
take the approach of using our photometry and existing K band photometry to compare the 
M stars among themselves to look for evidence of excess emission. Extrapolations from Ks 
to [24] require an accurate understanding of the photospheric colors that we determine in 
Section 13. 1[ We find no convincing evidence for any 24 /im excesses above the photospheric 
outputs. In section [X^ we use the ratio of our measured 70 /xm flux to a 70 /xm flux predicted 
from our 24 /xm measurement as our diagnostic tool. Again, we find no convincing evidence 
for a 70 /im excess at greater than the 3a level. The 160 /im upper limits are discussed in 
section 13.31 Finally, in sections 13.41 and 13.51 we compare the results on dust excesses around 
M stars with those obtained by previous Spitzer studies of main sequence stars in the solar 
neighborhood. 



3.1. Photospheric colors and 24 /im excess 

3.1.1. Photospheric colors 

Table [2] lists the Ks band magnitudes and the colors derived from t he infrared photona- 
etry. The band data were obtained mostly from the 2MASS catalog (jCutri. et al. 2003a ) 



but in some cases K measurements were obtained from other sources and transformed to 
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2MASS Ks (see notes to Table [2]). 

We show a color-color plot of K5 - [24] vs. V - K5 in Figure [H The V magnitudes 



were obtained from NStars (http://nstars.nau.edu/), where possible, with SIMBAD as a 



secondary source. Eleven of our 62 objects have no accurate V photometry available. The 
Ks - [24] colors are greater than zero, reflecting the fact that the spectral energy distributions 
of M stars do not transition to the Rayleigh- Jeans regime until well beyond 2 /im. They 
also become increasingly red with later spectral type, reflected by larger V - K5. With its 
long wavelength baseline, the V - color should be monotonic in spectral type/effective 
temperature. This plot should therefore allow us to identify excess candidate stars as outliers 
toward red K5 - [24] from the overall trend. In fact, there is one such object, V1581 Cyg B. 
However, the figure also shows significant scatter, substantially more than expected purely 
from the errors in V, K5 or [24], suggesting that there are significant variations in V - K5 
outside a purely monotonic relationship. 



3.1.2. Identification of Excesses Using IRFM for Stellar Temperatures 

Although accurate photometry for the entire sample is available in the near infrared, 
the issues shown in Figure 1 cannot be addressed with these measurements because M 
dwarf spectral type is not a sufficiently strong monotonic function of infrared colors. To 
i mprove over color-col or plots, we have used a modified form of the infrared fiux method 
(IBlackwell. et al.lll986l . and references therein) to determine effective temperatures for the 
stars in the sample, designated Tirfm- We can then plot K5 - [24] as a function of Tirfm 
rather than V - K5. 

The effective temperature can be expressed as: 



lOlog (Tirfm) = 42.32 - 5log 9.43 x 10^^ 



(24) 



2/3^ 



{JS,duY^\{T) (3{T)^ 



{niK + BCk) 



where a(T) is a term of order unity that is the ratio of the Planck function to the Rayleigh 
Jeans approximation at 24 /xm, and /5(T) is another term of order unity (/?(T) = 0.76 
independent of temperature for these stars) that is the ratio of the flux received through the 
24 fim fllter to the bolometric flux from the star, divided by the similar ratio for a blackbody 
at the effective t emperature of the s tar. F or stars of type M6 and later, we take bolometric 
corrections from iGolimowski. et al. J2004h . For stars earlier than M6, we used 



BCk = 0.42 + 1.45 {I - K) - 0.20 (/ - Kf 
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This expression is based on th at given bvlReid fc HawlevI (120001 ) . but has been modified to 
join better to the expression of iGohmowski. et aLl (120041 ) at the M5.5 to M6 transition. Our 
derived values of Tj^/m are hsted in Table [H 

Figure 2 shows the trend of K5 -[24] color vs. Tirfm- The scatter in the observed colors 
at a fixed temperature is small, consistent with it arising purely from the measurement 
errors. 

We have compared the colors of the low mass stars with synthetic colors genera ted from 



Brott fc Hauschildt 



20051 ) . These 



a grid of models to be published by Brott et al. (see 
models represent dusty atmospheres and are termed the "GAIA grid" because they were 
calculated in support of planning for the G AIA mission. They incorporate advances in 
the models published by lAUard. et al.l (I2OOII ). through improvements in parameters such as 
opacities and mixing lengths. The theoretical spectra were convolved with the wavelength- 
dependent instrumental response for the 2MASS K5 filter and the MIPS 24 //m filter. Similar 
convolutions were performed for a Kurucz model A-star spectrum, and the two were ratioed 
to obtain colors relative to A-star colors for the low mass stars. 

The synthetic colors are also compared with the measured ones in Figure [21 The obser- 
vations follow roughly the theoretical curve, although the observed colors are systematically 
up to 20% less red in the 2500 - 3500 K range. Among the contributors to this discrepancy 
could be the effects of dust in the atmospheres, which can modify the absorption feature 
strengths in the K-band as well as affect the radiation transfer and hence the emitted flux 
in the 24 /im band. These new observations should help provide a basis for improvement of 
models. 

Figure [2] also plots the positions of the stars using a conventional spectroscopic temper- 
ature calibration. We took spectral types from the NStars tabulation as of December 31, 
2006 wherever possible. In the remaining cases, the types are from SIMBAD as of the same 
date. In the latter instances, we consulted the primary references to verify the tabulated 
values. We used a type of M7.5 for LP 044-162, wher e the SIMBAD classiflcation of M6 5 
does not agree with the M7.5 found in the literature ( Gizis. et al. 200ol: Cruz, et al. 2003). 



From MO through M6, the temperature calibration is from lTokunagal (120001 ). while from M6.5 
through M9.5, we have assumed a linear relation with a slope of 200K per stellar sub-type. 
This relation joins the stellar and brown dwarf ter nperatures smoothly an d is consistent in 
approach with the temperature scales for L dwarfs (IBurgasser. et al.ll2002l ). Arrows indicate 
the shifts in effective temperature from these values to those determined through the IRFM 
in this work. If an object has a deflciency in flux at 24 fim relative to the conventionally 
determined temperature, then the IRFM calculation will indicate a higher temperature; if 
the object has a higher flux, then the IRFM will return a lower one. It can be seen that 
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most of the shifts are toward higher temperatures. 

If these stars had a combination of underestimated temperature plus an excess at 24 fim, 
then the corrected temperatures should show considerable scatter in the plot (because the 
bolometric corrections on which the plot depends are based on a broad range of wavelengths) . 
Instead, these points come to a well determined relationship with only small scatter. There- 
fore, we look for excesses among the stars where the IRFM shifts the temperature lower than 
the conventional calibration. An example is V1581 Cyg A , the star with a spectroscopic 
temperature of 2940K, adjusted to 2756K by the IRFM. However, this star has a low mass 
companion V 1581 Cyg C that is not resolved in our beam. This companion is faint enough 
in the vis ible not to bias the spectrurn significantly, but it can contribute significantly at 24 

(e.g., bawson fc de RobertisI boooh . A satisfactory fit to the properties of V 1581 Cyg 



A/C is a M6 primary and LI secondary. The predicted brightness difference between A and 
C is 3 magnitudes at R, while the predicted difference at [24] is only 0.8 magnitudes. The 
obse rved value is 3 .3 ma gnitudes at V (red version of V with an effective wavelength of 0.583 
iHenrv. et all Jl999h ): a difference of 0.8 magnitudes at [24] explains the behavior in the 
figure. Thus, the behavior in both bands is explained by the M6/L1 binary hypothesis. Al- 
though the excess in this case is not related to a debris disk, the temperature shift illustrates 
how a debris disk could be identified. 

The second object where we derive a lower temperature is V 1581 Cyg B, the single 
star suggested to have an excess in the color-color plot in Figure 1. It is not known to have 
a companion, but our results suggest it would be worthwhile to look for one. The third 
strongest example, LHS 3814, is also a double that is unresolved in our beam at 24 /im (we 
have plotted it at the temperature of the hotter member). The fourth strongest example, 
LP 944- 20, is thoug ht to be only ~300 Myr old, much younger than a typical object in the 



sample (jRibad 120031 ). Assuming it has lower surface gravity than the other members of the 
sample, the molecular bands in the 2 fim window will be strengthened consistent with its 
color deviation. All of the remaining objects fall in the region where no infrared excess is 
indicated. 



3.1.3. 24 l^m excesses 



Any star with a 24 fim excess should lie above the photospheric color locus in Figure 
2 or be shifted downward in temperature with our IRFM estimates. Although one could 
hypothesize that V1581 Cyg B has a debris-disk excess, since 24 fim excesses are very 
uncommon around FGK stars (iBryden. et al.ll2006l ) and since we find no large excess at 70 
Hm, this possibility is unlikely, and we will not consider it further. We therefore interpret 
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Figure [2] to show that there are no debris-disk-associated excesses at 24 /iin. 

We have repeated our analysis with artificial excesses at 24 /xm. An excess of 20% above 
the stellar photosphere results in a downward adjustment of the indicated stellar temperature 
by about 200 K. We can compare this value with the scatter in stellar colors. The thick solid 
line in Figure [2] can be fitted by 



Ks - [241 = -0.2915 + — — 

^ ^ 1.267 X 10-4 T+ (2.389 X 10-4 T)2 + (2.2 X 10-4 T)3 

A predicted 24 /im fiux density (P24 in Table [2]) was produced from this K5 - [24] and the 
measured K5. Then the ratio of the measured 24 /xm fiux density to P24 was formed for 
each star (F24/P24 in Table [2]). This ratio has a mean of 1.04 and a standard deviation of 
0.04 over our sample of 55 stars with derived Tj^/m- This standard deviation corresponds to 
about 100 K in temperature. Excepting the two cases discussed in the preceding section, all 
the stars in this sample either have spectroscopically determined temperatures lower than 
Tjr/m; or lu agreement with this temperature. Therefore, the limits on excesses are of the 
order of 20%, 2a. 



3.2. 70 /im Excesses 

Thirteen of the 20 stars with detected 70 /im fiuxes have high quality measurements 
(SNR > 4 and no interference from nearby objects or cirrus) that allow a good determination 
of 70 /im excess. We formed the ratio, given in Table [2] as F70/P70, of the measured 70 /im 
fiux to a predicted 70 /im fiux derived from a black body extrapolation of our measured 24 
/im fiux using Tirfm as the black body temperature. In Figure 3 we have plotted this ratio 
vs. Tirfm for these 13 stars. Objects with only photospheric emission appear at a value of one 
on this plot, while those with excess emission at 70 /im would appear at values above unity. 
None of the 13 stars shows a 70 /im excess larger than the 3a uncertainty in its (F70/P70) 
ratio. A 70 /im excess of 100% of the photospheric emission should be apparent in any of 
these stars. This null result is just at an interesting level of significance (see section 3.5.1 for 
discussion) . 

We also looked for 70 /im excesses by comparing our 70 /im detections among themselves. 
Any star with an excess should appear out of family with the, presumed, majority of non- 
excess stars. A ±3a range for F70/P70 of (0.59, 2.18) was obtained for the stars in Figure [3] 
by finding the mean and standard deviation of log{FiQ/Pjo). None of the 13 stars in Figure [3] 
exceed this range, confirming the conclusion of the previous paragraph. Among all twenty 70 
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/iin detections (SNR > 3.0) only AD Leo, YZ Get and LHS 39 fall outside of 0.59 < F70/P70 < 
2.18, the range indicating no significant excess. AD Leo is anomalously low and may indicate 
a problem with the measurement of this star. YZ Get may be contaminated by a background 
galaxy and should not be considered as having an excess without confirmation. LHS 39 is 
tantalizingly high, indicating the possibility of a 70 /im excess, but its 70 fj,m detection is 
only barely at 3.0 a significance and additional uncertainties entering the F70/P70 ratio result 
in too low a significance to claim an excess. 

3.3. 160 Upper Limits 

The 160 /zm band is more sensitive to the emission from colder dust, though the con- 
straints vary greatly from star to star depending on the quality of the 160 /im data. Lalande 
21185 and V1216 Sgr represent the range of 160 /xm limits within our data. We discuss 
the 160 /im upper limits along with the shorter wavelength measurements in the following 
section. 



3.4. Fractional dust luminosities 

Figure H] illustrates the limits on dust luminosity derived from our flux measurements. 
The upper limits on the ratio of dust luminosity to star luminosity (L^/L^,) are shown as 
a function of dust temperature for two stars, Lalande 21185 and V1216 Sgr, representing 
the range of limits obtained for stars detected at 70 /im. For each dust temperature, we 
calculate the maximum blackbody emission that is consistent with the observed 3a limits. 
Each of the curved segments in Figure H] corresponds to an individual observation at a specific 
wavelength, with the temperatures that an observation is most sensitive to depending on that 
wavelength. The 70 /xm observations are particularly useful for constraining the emission 
from dust with temperatures ~50 - 100 K, for which fractional luminosities greater than 
~10~^ are generally excluded. 

The constraints on L^/L^, in Figure H] can be translated into upper limits on dust mass, 
for an assumed dust size and orbital location. Assuming micron-size dust, the 160 /im upper 
limits correspond to dust mass limits of ~10~^ to 10"'' M^;^^^ of 20 K dust at 50 AU. 
Observations with moderate S/N at 70 /im translate to a maximum of ~10^'' MEarth of 50 
K dust at 5 AU, while 24 /im observations limit warm, 250 K dust to ~10-^ MEarth at 0.2 
AU (again assuming micron-size dust in both cases). 

While fractional dust luminosities similar to that of (3 Pic {Ld/L^ = 3 x 10^'^) would 
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have been detected in all stars with their photospheres detected at 70 /im, luminosities as 
low as that of Vega [La/L^,^ 2 x 10'^) could only be detected around stars with optimal 
data such as Lalande 21185. 



3.5. Comparison of M stars with F, G, and K Stars 



3.5.1. Previous Work 



Debris disk excesses in M stars have been elusive, partly because of the increase in 
dete ction difficulty w i th de creasing stellar luminosity (see Section [3.5.2p . ISong. et al.l (120021 ) 
and iPlavchan. et al. J2005h have searched for excesses from M stars from IRAS data com- 
plemented by ground based measurements. Both studies find no credible excesses around 
mature M stars, that is those not associated with known young clusters nor showing spec- 
troscopic signs of youth. Our results confirm this con clusion, using the more sensitive and 
accurate Spitzer data. Using reprocessed IRAS data, iBackman. et al.l (119861 ) . reported an 
infrared excess from Ross 128 at 60 and 100 /im. This was not confirmed by subsequent 
processing in the IRAS Faint Source Survey, and is also contradicted by the Spitzer/MlPS 
results. The 0.12 Jy excess they reported at 60 fj,m is 8 times larger than our 3 sigma upper 
limit at 70 um, an d is therefore confidently excluded by our data. Using Spitzer spectroscopy. 
Jura, et al.l (120041 ) report no excess for Ross 128. This leaves no known example of a mature 
main sequence M dwarf known to possess a far-infrared excess. 

It is interesting to compare the present results with similar studies of stars of other 
spectral types to see whether M star debris disks are less common than disks around earlier 
type stars. As a comp arison sample, we ha ve used both the study of F5-K5 ma in-sequence 



f ield s tars discussed by lBryden. et al.l (120061 ) and the similar set of stars studied by lKim. et al. 



( I2OO5I ). We combine all the data without regard to whether they apply to single stars, 
binaries, or stars with radial velocity planets. There are 69 stars in Bryden et al. and 35 in 
Kim et al., for a total of 104. Of these stars, 12 have detected excesses. We need to add ~2 
more detections to allow for the weighted con tribution of the four stars with bright excesses 
detected by IRAS and not re-observed (see iBryden. et al.l 120061 ) . We find that 13 ± 4% 



of the se stars have excesses greater than 100%. Only one of the 69 stars in iBryden. et al. 



( I2OO6I ). HD 69830, shows evi dence for a signifi cant 24 fim photometric excess more than 
15% above the photosphere. iKim. et al.l (120051 ) included only stars with cold excesses, so 
we cannot derive an independent statistic in this regard. However, it is clear that 24 /im 
excesses are uncommon in solar type main sequence field stars. We find no convincing 24/im 
excess in a sample of 56 M stars. 
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The current sample of 20 M dwarfs with 3cr or better 70 fim detections is quite small 
for purposes of drawing confident statistical inferences of debris disk frequency. Binomial 
statistics indicate that our result of no disk detections in the 20 stars only allows us to derive 
an upper limit of 14%, with 2a confidence, to the incidence of disks among M dwarfs within 
5pc of the Sun. This upper limit is right at the measured incidence level f or more massive 
stars. Our upper limit is also consistent with a value of 13% reported by iLestrade. et al. 
( I2OO6I ) for 20 to 200 Myr old field M dwarfs. However, we note that our 5 pc sample is likely 
to contain a large fraction of stars older than 200 Myr and may not be directly comparable. 



Song, et al.l (|2002| ) find that a number of extremely young M stars haye det ectable 



excesses (i.e.. Hen 3-600 and AU Mic, both less than 15 M yr old). ILow. et al.l (120051 ) report 
sensitive Spitzer measurements in the young TW Hya association wh ere, of the 2 4 stars 
observed, 5 have significant excess at 24 /im and 6 have excess at 70 /im. I Chen, et al.l (120051 ) 
confirm the excess in AU Mic. Thus, excesses are not uncommon in late-type stars, including 
M-type ones, at an age of ~10 Myr. However, they virtually completely disappear in mature 
low-mass stars. 



3.5.2. Approach for Quantitative Comparison 

To expand on the results in the preceding section, we need to use an appropriate metric 
to characterize debris systems. We will adopt the 70 /im fractional excess, ^ excess photosphere^ 
where ¥ excess is the flux density of the infrared excess and Fphotosphere that of the star itself, 
both at 70 /im. We use the 70 /xm fractional excess as a proxy for the fractional luminosity, 
Liji/Lgtar, because it should be representative of disk luminosity for objects with no excess 
at shorter wavelengths. 

As stars become less luminous and of lower temperature, it becomes increasingly difficult 
to measure infrared excesses to a given limit in fractional luminosity. This issue can be 
understood by considering a set of stars with varying temperature. The total luminosity of 
these objects, and thus the dust luminosity, will vary as T'^R^ (where R is the stellar radius), 
but their photospheric outputs in the far infrared will be in the Rayleigh- Jeans regime and 
will only go as TR^. Thus, for a given threshold for excess detection above the photosphere, 
the fractional luminosity threshold will go as T~^. 
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3.5.3. Application to M Star Sample 

To obtain a quantitative measure of the average level of excess in the M star sample, we 
proceed in a manner analogous to stacking, a technique often used to enhance the average 
signal to noise ratio with observations of a suite of similar sources. Our various M stars 
have different expected signal levels and, as a result, the measurements are of intrinsically 
differing ratios of signal to noise. Therefore, we stack by taking the weighted average of the 
measured excess ratios at 70/im. We eliminate from this average all objects with a nominal 
SNR < 4. We adopt for the remaining sources a net error consisting of the statistical error 
indicated in Table 1, plus a systematic erro r of 6% to reflect th e scatter in measurements 



of bright objects at 70 /im with MIPS (e.g. iGordon. et al.ll2005l ). The result is a weighted 



average excess ratio at 70 fim of 1.02 ±0.05; that is, on average, there is no excess from these 
stars. We place a 2a upper limit of 1.12 to the average excess ratio, or 0.12 to the excess 
(relative to the photospheric emission). We also compute an average effective temperature 
for this sample as a weighted average of the individual stellar effective temperatures where 
the weights are those used for the excess ratios. The result is T = 3555 K. 



As a comparison sample, we use the FGK stars in lBryden. et al.l (120061 ). We add in the 
four additional stars listed by Bryden et al. but measured in other Spitzer programs at half 
weight, since Bryden et al. are reporting on about half of their sample. We also add in DY 
Eri and eta Gas from this program. The result is an average excess ratio of 1.37 ± 0.033. 
The average temperature of this sample computed in the same way as for the M stars is 5763 
K. If we scale the average excess for this sample as T^, we predict that the average for the 
M stars if they have the same fractional luminosity in infrared excess emission is 0.09, just 
within the upper limit from our observations of 0.12. 

Another effect that might lead to more prominent debris excesses in M stars is the scaling 
of dust mass to produce a given excess level, as a function of stellar luminosity. For a given 
fractional luminosity upper limit at a given wavelength, the minimum mass limits derived 
from the M star observations reported here are lower than would be the corresponding limits 
from more luminous stars. The reason is that the fractional dust luminosity is proportional to 
the total solid angle subtended by the dust as viewed from the star and that the temperature 
of the dust increases with the star's luminosity but decreases with distance from the star. 
Thus, at a fixed dust temperature, an increase in stellar luminosity requires that the dust be 
more distant which, in turn, requires more dust to subtend enough solid angle to keep L^/L^, 
constant. In the simple case of a single fixed grain size, the number of grains to achieve a 
certain fractional luminosity at a given temperature scales with the luminosity of the star. 
Thus the minimum dust mass associated with an M star with ~ 0.01 Lq solar luminosities 
would be ~ 1% of the corresponding minimum mass associated with a solar mass star, again 
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assuming the same limit to the fractional luminosity at a given wavelength. For typical 
M stars in the present sample, therefore, it is not surprising to note that the limits on 24 
/im excesses correspond to masses as low as 10~^ earth masses, while at 70 yum, the non- 
detections could imply masses b elow 10"'' earth mas ses. Similar Spitzer/MIPS surveys for 
IR excess around FGK stars (e.g. iBryden. et al.ll2006l ) can detect dust with lower luminosity 
relative to the central star (Ldust/^'star as low as 5 x 10~^), but are less sensitive in terms of 
overall dust mass. The corresponding limit for non-detections of 70 /im excess around FGK 
stars is > 10^^ ^Earth, an order of magnitude above that for M stars. 



3.5.4- Discussion 

The fractional luminosity is the appropriate metric if we assume that hypothetical debris 
disks would have the same covering fraction independent of stellar type, thus absorbing and 
reradiating in the far infrared an equivalent fraction of the stellar output. Although the 
fractional luminosity provides a reasonable first-order metric to compare excesses in different 
stellar populations, there are many other considerations. For example, one would expect that 
cool, low-mass stars might be less effective at heating circumstellar dust than hotter ones, 
even after allowing for the difference in luminosity. The degree to which this expectation is 
realized depends on the optical properties of the debris disk particles. They are likely to be 
relatively large grain aggregates, whose behavior may differ substantially from that of typical 
interstellar grains. The grain loss mechanisms - photon pressure, Poyting-Robertson drag, 
and corpuscular winds - are generally different in low mass stars than in high mass ones, and 
may cause significant differences in their disks (e.g., strong winds may destroy them quickly). 
Also, our study probes debris disks at 70 fim and hence at a temperature of about 50 K. 
The equilibrium distance from a star of luminosity L for material at a given temperature 
scales as the square root of L. Given the steep dependence of stellar luminosity on mass, we 
are probing distinctly different physical regions around the M stars in this sample than the 
regions probed in observations of solar-like stars. 



4. Conclusions 

We observed a sample of 62 M stars to search for debris disk dust and to study photo- 
spheric models. We demonstrate a modified infrared flux method to determine the effective 
temperatures of these stars, taking advantage of the high accuracy measurements available 
at 24 fim with Spitzer. Using these effective temperatures, the K5 - [24] colors of the stars 
agree reasonably well with the theoretical models of Brott et al. (private communication). 
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None of the stars in our sample has a high- weight detection of a mid- or far-infrared excess. 
In particular, the sample lacks the occasional stars with large excess that typifies samples 
of more massive stars. We compute a weighted average of the observed excess ratios at 70 
fim. as a metric to characterize the excesses in this and other samples of stars. We find that 
the average excess in the M stars is at least four times less than the average in a sample 
of FGK stars. This limit is just at the level for equal fractional infrared luminosity in the 
two samples. We have set low hmits (10~^ to 10~^ ^Earth depending on location) for the 
maximum mass of dust possible around our stars. 

This work is based, in part, on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, Cahfornia Institute of Technology un- 
der contract with NASA. Support for this work was provided by NASA through contract 
1255094. This work also makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center /California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. We would further like to thank 
Peter Hauschildt for supplying us with results from the latest GAIA model spectra of late 
type stars. 
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Table 1. Observational Data 



Name 


Gliese 


LHS 


SpTi 


F24(mJy) 


F7o(mJy)* 


Fi6o(mJy,; 


AX Mic 


GJ 825 


66 


MO 


504.00 ± 0.57 


61.8 ± 3.8 


<75 


LHS 38 


GJ 412 A 


38 


Ml 


123.94 ± 0.21 


15.8 ± 3.9^ 


<67 


GX And 


GJ 15 A 


3 


Ml. 5 


249.45 ± 0.29 


27.2 ± 2.2 


<102 


Kaptcyn's Star 


GJ 191 


29 


Ml. 5 


98.46 ± 0.14 


4.6 ± 2.4*^ 




LHS 70 


GJ 887 


70 


Ml. 5 


436.96 ± 0.43 


47.3 ± 3.9 


<53 


Lalande 21185 


GJ 411 


37 


M2 


481.55 ± 0.40 


43.0 ± 3.5 


<62 


LHS 1 


GJ 1 


1 


M3 


159.13 ± 0.28 


16.2 ± 3.2 




AD Leo 


GJ 388 


5167 


M3 


151.62 ± 0.20 


8.4 ± 2.3 




V2306 Oph 


GJ 628 


419 


M3 


104.83 ± 0.23 


19.2 ± 2.8 


<154 


T TTO A Art 

LHb 449 


GJ 674 


449 


M3 


120.51 ± 0.22 


10.4 ± 7.3 


<963 


LHb 450 


GJ 687 


450 


M3 


160.17 ± 0.21 


28.1 ± 3.4 


<83 


LHS 58 


GJ 725 A 


58 


M3 


183 ± 3 


31.9 ± 3.3"^ 


<73'^ 


LHS 3685 


GJ 832 


3685 


M3 


163.15 ± 0.18 


19.5 ± 4.1 


<57 


T TTO m -1 A 

LHS 3814 


T a^t\ A 

GJ 860 A 


3814 


M3 


t A ^ f\ 1 A 

146.0 ± 0.4 


117 ± 61" 




GQ And 


GJ 15 B 




M3.5 


49.46 ± 0.20 


5.3 ± 1.9^ 




Luyten's Star 


GJ 273 


33 


M3.5 


127.16 ± 0.22 


11.4 ± 3.7 


<71 


LHS 59 


GJ 725 B 


59 


M3.5 


115 ± 3 


11.0 ± 3.3"* 




Vlilb bgr 


GJ 729 


3414 


M3.5 


83.67 ± 0.27 


17.0 ± 5.7s 


<318 


IL Aqr 


GJ 876 A 


530 


M3.5 


108.27 ± 0.16 


9.2 it 2.2 


<67 


Ross 128 


GJ 447 


315 


M4 


66.32 ± 0.24 


7.5 ± 5.0 


<61 


Barnard's Star 


GJ 699 


57 


M4 


188.45 ± 0.25 


29.1 ± 4.4 


<141 


YZ Get 


GJ 54.1 


138 


M4.5 


33.69 ± 0.19 


9.0 ± 2.0 




TZ Ari 


GJ 83.1 


11 


M4.5 


27.71 ± 0.21 






DY Eri C 


GJ 166 G 


25 


M4.5 


49.05 ± 0.14 


4.0 ± 2.0 




T rcr'^T TV T A 

V577 Mon A 


/~i T A A 

GJ 234 A 


1849 


M4.5 


77.97 ± 0.20 


12.2 ± 1.9' 


<128 


LHS 451 


GJ 682 


451 


M4.5 


62.15 ± 0.23 


k 




EZ Aqr 


GJ 866 A 


68 


M5 


78.73 ± 0.19 


-4.2 ± 3.5 




LHS 39 


GJ 412 B 


39 


M5.5 


9.89 ± 0.16 


10.1 ± 3.4' 




FL Vir 


GJ 473A 


333 


M5.5 


49.4 ± 0.4 


3.0 ± 2.0 




Proxima Cen 


GJ 551 


49 


M5.5 


235.06 ± 0.37 






Ross 248 


GJ 905 


549 


M5.5 


52.78 ± 0.16 


10.7 ± 2.0° 




LHS 2 


GJ 1002 


2 


M5.5 


14.60 ± 0.18 


4.1 ± 2.6 




LHS 1565 


GJ 1061 


1565 


M5.5 


29.65 ± 0.22 


p 


<63 


V1581 Cyg A 


GJ 1245 A 


3494 


M5.5 


27.0 ± 1.6 


-12 ± 10^ 


<387 


LHS 1375 


GJ 3146 


1375 


M5.5 s 


3.31 ± 0.13 






LHS 288 


GJ 3618 


288 


M5.5 


10.71 ± 0.14 


r 




UV/BL Get 


GJ 65 B 


10 


M6 


95.74 ± 0.22' 


5.6 ± 1.9= 


<68 


GL 316.1 


GL 316.1 


2034 


M6 s 


1.35 ± 0.11 






Wolf 359 


GJ 406 


36 


M6 


49.90 ± 0.30 


t 




V1581 Cyg B 


GJ 1245 B 


3495 


M6 


18.0 ± 0.7 


-12 ± 101 




GJ 283B 


GJ 283 B 


234 


M6.5 s 


1.05 ± 0.13 






DX Cnc 


GJ 1111 


248 


M6.5 


17.18 ± 0.22 


6.0 ± 3.0" 


<101 


LHS 292 


GJ 3622 


292 


M6.5 


9.41 ± 0.13 


1.1 ± 1.9^ 




LHS 2930 


GJ 3855 


2930 


M6.5 s 


1.63 ± 0.1 






GJ 644C 


GJ 644 G 


429 


M7 


4.17 ± 0.12 






LHS 3003 


GJ 3877 


3003 


M7 


3.75 ± 0.11 
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Table 1 — Continued 



Name 


Gliese 


LHS 


SpT^ 


F24(mjy) 


F7o(mJy)* F160 (niJy,3(T) 


2MASS J09522188-1924319 






M7s 


0.53 ± 0.07 




2MASS J23062928-0502285 






M7.5 s 


1.35 ± 0.12 




LP 044-162 






M7.52 


1.08 it O.OS 




LP 349-025 






MS s 


2.31 ± 0.12 




LP 326-021 






MS s 


0.S9 it 0.09 




LP 412-031 






MS s 


0.74 it 0.07 




LP 771-021 






MS s 


0.41 it 0.05 




LHS 2397 A 


GJ 3655 


2397A 


MS s 


0.99 it 0.06 




ESO 207-61 






MS s 


0.16 it 0.04 




TVLM 513-46546 






MS.5 s 


1.03 it 0.1 




LP 944-20 






M9 s 


3.95 it 0.05 


7 ± 3W 


LHS 2065 


GJ 3517 


2065 


M9 


1.9 it 0.09 




TVLM 868-110639 






M9 s 


0.68 it 0.03 




BRI 0021-0214 






M9.5 s 


1.26 ± 0.05 




2MASSW J1733189-I-463359 






M9.5 s 


0.28 it 0.03 




DENIS-P J0021.0-4244 






M9.5 s 


0.26 it 0.04 





* Entries with an ellipsis ( • • ■ ) in the 70 fim column are from Spitzer program 56. 
Spectral types are from the NStars data base except for those marked "s" which are from Simbad. 
^See text for LP 044^162 spectral type source. 
''LHS 38: Some cirrus at 70 /im 

''Kapteyn's: No source apparent at 70 /im. Photometry at reference position. 
"^LHS 449: Bad cirrus confusion at 70 ^m 

"^LHS 58 & 59: 14" separation. 70 ^m fluxes from deblending. 160 //m upper limit is for the combination. 
''LHS 3814: Bad cirrus, especially at 70 ixm 

'GQ And: 35" separation from GX And: possible slight blending at 70 /um. 
SV1216 Sgr: Bad Cirrus at 70 fira 

^YZ Get: Diffuse source, probable background contamination 
'TZ Ari: Bright background source 20" north a<;counts for flux 

jV577 Mon A: 2" separation from GJ23B. Photometry includes both sources in both bands 

''LHS 451: Bad cirrus confusion: no source apparent at 70 /im 

'LHS 39: 70 /im photometry done at 24 /im position. 32" separation from LHS 38 
™Proxima Gen: Very bad cirrus at 70 /im. 

"Ross 248: 70 /tm source offset slightly from 24 /im source; possible contamination. 
PLHS 1565: Bright background galaxy. No 70 /tm photometry possible 

1V1581 Gyg A & B: S" separation and bad cirrus. 70 /im photometry includes both sources 

"^LHS 288: Not on the 70 /tm image (error in requested pointing) 

''UV Cet/BL Get: 2" separation. Photometry is for both sources in both bands 
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*Wolf 359: Not on the 70 /im image (error in requested pointing) 

"DX Cnc: No source apparent at 70 ^m. Photometry at reference position. 
^LHS 292: 70 um photometry at 24 ^m position. Nearby 70 ^m source. 
"LP 944-20: Faint source at 70 ^J,m■, another 70 /im source 50" north 



Table 2. Colors and Derived Quantities 



Name 


Dist(pc) 


SpT 






B.C. 


Ks* 


Ks-[24]t 


P24* 


F24/P24* 


P70* 


F70/P70* 


AX Mic 


3.95 


MO 


3800 


3987 


2.370 


3.07 ±0.02^ 


0.19 ± 0.06 


502.04 


1.00 


58.37 


1.06 


LHS 38 


4.85 


Ml 


3680 


3591 


2.530 


4.77 ±0.02 


0.37 ± 0.06 


116.58 


1.06 


14.44 


1.09 


GX And 


3.56 


Ml. 5 


3605 


3646 


2.470 


4.02 ± 0.02 


0.38 ± 0.06 


228.78 


1.09 


29.03 


0.94 


Kapteyn's Star 


3.92 


Ml. 5 


3605 


3554 


2.474 


5.05 ± 0.02 


0.40 ± 0.06 


88.49 


1.11<* 


11.46 


< 0.6"* 


LHS 70 


3.29 


Ml. 5 


3605 


3636 


2.530 


3.36 ± 0.02* 


0.33 ± 0.06 


421.41 


1.04 


50.86 


0.93 


Lalande 21185 


2.54 


M2 


3530 


3469 


2.530 


3.35 ± 0.02*' 


0.42 ± 0.06 


432.61 


l.ll'^ 


56.10 


0.77'' 


LHS 1 


4.36 


M3 


3380 


3563 


2.530 


4.52 ± 0.2 


0.39 ±0.21 


148.05 


1.07 


18.55 


0.87 


AD Leo 


4.89 


M3 


3380 


3360 


2.648 


4.59 ± 0.02 


0.41 ± 0.06 


148.76 


1.02 


17.73 


0.47 


V2306 Oph 


4.24 


M3 


3380 


3265 


2.710 


5.08 ± 0.02 


0.49 ± 0.06 


99.15 


1.06 


12.28 


1.56 


LHS 449 


4.54 


M3 


3380 


3471 


2.580 


4.86 ± 0.02 


0.42 ± 0.06 


112.58 


1.07 


14.07 


< 1.6 


LHS 450 


4.54 


M3 


3380 


3332 


2.710 


4.55 ± 0.02 


0.43 ± 0.06 


155.96 


1.03 


18.74 


1.50 


LHS 58 


3.53 


M3 


3380 


3369 


2.650 


4.43 ± 0.02 


0.45 ± 0.06 


171.81 


1.07 


21.40 


1.49 


LHS 3685 


4.93 


M3 


3380 


3657 


2.470 


4.47 ± 0.02^= 


0.37 ± 0.06 


150.66 


1.08 


18.98 


1.03 


LHS 3814 


4.03 


M3 


3380 


3211 


2.648 


4.78 ± 0.03 


0.56 ± 0.06 


132.40 


1.10 


17.13 


< 10.7 


GQ And 


3.56 


M3.5 


3280 


3157 


2.708 


5.95 ± 0.02 


0.55 ± 0.06 


46.12 


1.07 


5.81 


< 1.0 


Luyten's Star 


3.79 


M3.5 


3280 


3273 


2.710 


4.86 ± 0.02 


0.49 ±0.06 


119.93 


1.06 


14.90 


0.77 


LHS 59 


3.53 


M3.5 


3280 


3239 


2.710 


5.00 ± 0.02 


0.52 ± 0.06 


106.87 


1.08 


13.48 


0.82 


V1216 Sgr 


2.97 


M3.5 


3280 


3215 


2.710 


5.37 ± 0.02 


0.54 ± 0.06 


76.76 


1.09 


9.81 


< 1.8 


IL Aqr 


4.70 


M3.5 


3280 


3240 


2.760 


5.01 ± 0.02 


0.46 ± 0.06 


105.85 


1.02 


12.69 


0.72 


Ross 128 


3.35 


M4 


3180 


3130 


2.760 


5.65 ± 0.02 


0.57 ± 0.06 


61.53 


1.08 


7.79 


< 1.9 


Barnard's Star 


1.83 


M4 


3180 


3039 


2.850 


4.52 ± 0.02 


0.57 ± 0.06 


181.72 


1.04 


22.19 


1.31 


YZ Get 


3.72 


M4.5 


3105 


3052 


2.810 


6.42 ± 0.02 


0.60 ± 0.06 


31.38 


1.07 


3.97 


2.27 


TZ Ari 


4.45 


M4.5 


3105 


2995 


2.810 


6.65 ± 0.02 


0.62 ± 0.06 


26.11 


1.06 


3.27 




DY Eri C 


5.03 


M4.5 


3105 


3058 


2.810 


5.96 ± 0.02 


0.55 ± 0.06 


47.80 


1.03 


5.77 


< 1.0 


V577 Mon A 


4.09 


M4.5 


3105 


3074 


2.810 


5.49 ± 0.02 


0.59 ± 0.06 


73.14 


1.07 


9.17 


1.33 


LHS 451 


5.01 


M4.5 


3105 


3302 


2.710 


5.61 ± 0.02 


0.46 ± 0.06 


59.43 


1.05 


7.28 




EZ Aqr 


3.45 


M5 


3030 


2855 


2.943 


5.54 ±0.02 


0.65 ± 0.06 


78.21 


1.01 


9.31 


< 1.1 


LHS 39 


4.85 


M5.5 


2940 


2811 


2.943 


7.84 ± 0.03 


0.69 ± 0.06 


9.65 


1.03 


1.17 


< 9.0 


FL Vir 


4.39 


M5.5 


2940 


2866 


2.943 


6.04 ± 0.02 


0.63 ± 0.06 


49.04 


1.00 


5.80 


< 1.0 


Proxima Gen 


1.30 


M5.5 


2940 


2855 


2.943 


4.35 ± 0.02^ 


0.64 ± 0.06 


234.04 


1.00 


27.80 




Ross 248 


3.16 


M5.5 


2940 


2928 


2.940 


5.93 ± 0.02 


0.60 ± 0.06 


52.46 


1.01 


6.23 


1.72 


LHS 2 


4.69 


M5.5 


2940 


2792 


2.943 


7.44 ± 0.02 


0.72 ± 0.06 


14.10 


1.04 


1.73 


< 4.5 


LHS 1565 


3.69 


M5.5 


2940 


2879 


2.940 


6.61 ± 0.02 


0.66 ± 0.06 


28.80 


1.03 


3.50 




V1581 Gyg A 


4.54 


M5.5 


2940 


2756 


2.940 


6.85 ± 0.02 


0.79 ± 0.08 


24.82 


1.09 


3.20 


< 10.3 


LHS 1375 


8.50 


M5.5 


2940 


2821 


3.034 


8.98 ± 0.05 


0.65 ± 0.09 


3.36 


0.99 







Table 2 — Continued 



NamG 




op 1 




T I- 
^ irjm 


B.C. 






-T 24 


-T 24/-r 24 


r 70 


r 70/r70^ 


LHS 288 


4.49 


M5.5 


2940 


2832 


2, 


.943 


7.73 ± 0.03 


0.67 ± 0.06 


10.54 


1.02 


1.27 




UV/BL Get 


2.68 


M6 


2850 


2874 


2, 


.940 


5.34 ± 0.02 


0.66 ± 0.06 


93.03 


1.03 


11.32 


< 0.5 


GL 316.1 


14.06 


M6 


2850 


2724 


3, 


.034 


10.05 ± 0.05 


0.74 ± 0.12 


1.33 


1.02 






Wolf 359 


2.39 


M6 


2850 


2857 


2, 


.943 


6.08 ±0.02 


0.69 ± 0.06 


47.51 


1.05 






V1581 Gys B 


4.54 


M6 


2850 


2602 


3, 


.030 


7.39 ±0.02 


0.88 ± 0.07 


16.89 


1.07 


2.14 


< 15.4 


GJ 283B 


8.90 


M6.5 


2750 


C 






9.26 ±0.05 


-0.32 ± 0.15 


2.48 


0.42<= 






DX Gnc 


3.63 


M6.5 


2750 


2746 


3, 


.030 


7.26 ±0.02 


0.71 ±0.06 


17.12 


1.00 


2.04 


< 4.4 


LHS 292 


4.54 


M6.5 


2750 


2675 


3, 


.069 


7.93 ± 0.03 


0.73 ± 0.06 


9.67 


0.97 


1.12 


< 5.1 


LHS 2930 


9.63 


M6.5 


2750 


2711 


3, 


.069 


9.79 ± 0.05 


0.69 ± 0.10 


1.70 


0.96 






GJ 644G 


6.45 


M7 


2650 


2700 


3, 


.060 


8.82 ±0.05 


0.74 ±0.08 


4.19 


0.99 






LHS 3003 


6.40 


M7 


2650 


2704 


3, 


.060 


8.93 ± 0.05 


0.73 ± 0.08 


3.78 


0.99 






2MASS J0952-1924 




M7 


2650 


C 






10.87 ±0.05 


0.55 ±0.16 


0.67 


0.79"= 






2MASS J2306-0502 




M7.5 


2550 


2479 


3, 


.070 


10.30 ± 0.05 


0.99 ± 0.12 


1.26 


1.07 






LP 044-162 


12.5 


M7.5 


2750 


2589 


3, 


.070 


10.40 ± 0.05 


0.85 ± 0.11 


1.06 


1.02 






LP 349-025 


7.8 


M8 


2450 


2570 


3, 


.100 


9.57 ±0.05 


0.84 ± 0.09 


2.30 


1.00 






LP 326-021 




M8 


2450 


2559 


3, 


.100 


10.62 ± 0.05 


0.86 ±0.13 


0.88 


1.01 






LP 412-031 


14.64 


M8 


2450 


2728 


3, 


.056 


10.64 ±0.05 


0.68 ±0.13 


0.77 


0.96 






LP 771-021 


16.23 


M8 


2450 


2591 


3, 


.100 


11.42 ± 0.05 


0.82 ±0.15 


0.41 


0.99 






LHS 2397 A 


14.29 


M8 


2450 


2380 


3, 


.100 


10.74 ±0.05 


1.09 ±0.10 


0.92 


1.08 






ESQ 207-61 


18.48 


M8 


2450 


C 






12.11 ± 0.05 


0.49 ±0.28 


0.25 


0.65« 






TVLM 513-46546 


9.82 


M8.5 


2350 


2404 


3, 


.120 


10.71 ± 0.05 


1.11 ±0.13 


0.92 


1.12 






LP 944-20 


4.97 


M9 


2250 


2123 


3, 


.150 


9.55 ±0.02 


1.41 ±0.06 


3.58 


1.10 


0.48 


< 18.8 


LHS 2065 


8.53 


M9 


2250 


2420 


3, 


.150 


9.94 ±0.05 


1.00 ±0.09 


1.85 


1.03 






TVLM 868-110639 


17.39 


M9 


2250 


2193 


3, 


.150 


11.35 ± 0.05 


1.30 ±0.09 


0.63 


1.08 






BRI 0021-0214 


12.12 


M9.5 


2150 


2284 


3, 


.180 


10.54 ±0.05 


1.16 ±0.08 


1.21 


1.04 






2MASSW J1733+4633 




M9.5 


2150 


2512 


3, 


.180 


11.89 ± 0.05 


0.87 ±0.14 


0.28 


0.99 






DENIS-P J0021. 0-4244 




M9.5 


2150 


d 






12.30 ±0.05 


1.20 ±0.18 


0.22 


1.21'^ 







to 
00 



*Kfi is the 2MASS magnitude. K magnitudes footnoted and were converted to Kg using the 2MASS-GIT transformation in lGutri. et al.l 
ll2003bl) 

t24um zero point magnitude taken as 7.14 Jy. 

•'■Predicted values of the flux at 24 and 70 fim, P24 and P70 are derived as described in the text, except that footnoted values use Teff instead 
of Tirfm Upper limits are at the 2 cr level. 

'^Photometry from Mould, J.R. & Hyland, A.R., 1976, ApJ, 308, 399 



''Photometry from Persson et al., 1977, A J, 82, 9 

'^Photometry from Lcggett, S.K. 1992, ApJS, 82, 351 

'^Excluded from K5-[24] calculation due to low SNR, but fits the trend. 

"^Excluded from Kg-[24] calculation due to low SNR, and is highly discrepant. 
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Fig. 1. — K5-[24] colors vs. V - K5. The large symbol is for V1581 Cyg B, the only object 
with a possible excess at 24//m identified through this test. 
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Fig. 2. — K5'-[24] colors vs. Teff and Tirfm- The open circles show Teff from the spectral 
type, while the filled boxes are the result of deriving Tirfm from the infrared fiux method, as 
described in the text. Unresolved binaries are indicated by larger circles and young objects by 
open squares. Cases with large unexplained differences are connected with arrows from the 
spectral estimate to the infrare d flux method one. The th in solid line shows the prediction 
from the atmospheric models of iBrott fc Hauschildtl (120051 ). The thick solid line is the flt to 
K5-[24]vs. Tirfm derived in this paper. The dashed line shows the color of a black body. 
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Fig. 3. — The ratio of measured to predicted 70/im flux vs. the effective temperature of the 
sources determined from the infrared flux method. Only sources that were unambiguously 
detected at 70 fim are included (i.e. SNR > 3 and no obvious indications of contamination 
from background sources; see notes to Table [2]). The 70 fim prediction is extrapolated from 
the measured 24 fim flux using a black body at Tj^/m for each star, as described in the text. 
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Fig. 4. — Detection limits for circumstellar dust. The MIPS detection levels are shown for 
two stars as a function of dust luminosity and temperature. The light grey region is ruled 
out by the detection upper limits for Lalande 21185 (solid line), a relatively bright star with 
a clean 160 fim image, while the dark grey region is ruled out for V1216 Sgr (limits shown as 
a dotted line), a star located in a very noisy 160 /^m field. Both stars have la errors of ~5% 
at 24 iim, but their relative accuracies diverge at longer wavelengths; la errors at 160 /xm 
are 0.9 and 70 times the stellar photosphere for Lalande 21185 and V1216 Sgr respectively. 
In both cases, 3 a limits are plotted. 



